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Abstract 

We generalize the exactly solvable corner growth models by choosing the rate of the exponential 
distribution at + hj and the parameter of the geometric distribution aibj at site where 

and (bj)j^i are jointly ergodic random sequences. We identify the shape function in terms of a 
simple variational problem, which can be solved explicitly in some special cases. 


1 Introduction 


The corner growth model is a freqnently studied model of random growth. See [29j for a detailed 
introduction, and [22], [28] for an overview of related research. The model describes a cluster of 
sites that emerges from the corner of a quadrant and grows over time. We represent the quadrant 
with and the cluster with a family of subsets St cz indexed by time t ^ 0. Each site 
(i, j) e waits for {i — l,j) if i > 1 and {i,j — 1) if j > 1 to be in the cluster and, after an 
additional waiting time of W{i,j), joins the cluster permanently at time G{i,j). More precisely, 
St = {(bi) £ ^ t} for t 5= 0 and 

Gii,j) = G{i-l,j)vG{iJ-l) + W{i,j) fori,jeN, (1.1) 

where G{i, 0) = 0(0, j) = 0 for i,j e N. As first observed in [53], the preceding recursion implies 
G{m,n) = max V W(i,j) form,nGN, (1.2) 

ttpEIi 1 ™ 


where is the set of all directed paths from (m, v) to (u', u'), that is, all finite sequences tt = 

((4, jfe))issfcss/ inZ^ such that {ik+i-ik, jk+i-jk) e {(1,0), (0,1)} for 1 ^ fc < Z, = {u,v) and 

(hGi) = W ,v'). We will refer to the random quantities {W{i,j) : b j e N} and {G{i,j) : b j e N} as 
weights and last-passage times, respectively, because of connection with directed last-passage 
percolation. The problem is typically to understand the statistical properties of the last-passage 
times given the joint probability distribution P of the weights. 

Consider the corner growth model in which the weights are i.i.d. and exponentially distributed 
with rate A > 0 i.e. P{W{i,j) ^ x) = for i,j e N and x ^ 0. We will refer to this special 
case as the exponential model. In a seminal paper [25], H. Rost observed the equivalence of the 
exponential model with the totally asymmetric simple exclusion process (TASEP) and proved that 


lim 
n—*-00 


G([nsJ,[ntJ) _ (vs-f Vt)^ 


A 


for s,t > 0 P-a.s. 


(1.3) 


He also interpreted (1.3) as an asymptotic shape result in the sense that the rescaled cluster Sn/n 
converges as n ^ oo to the parabolic region {s,t e IR+ : ^/s -)- s/t V VA} in the Hausdorff metric. 
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The discrete counterpart of the exponential model is the geometric model in which the weights are 
i.i.d. and geometrically distributed with (fail) parameter q e (0,1) i.e. P{W{i,j) ^ k) = q'^ for 
i,j G N and k G For the geometric model, one can also compute that 


lim 

n—»-GO 


G{[ns\,[nt\) _ q 


{s + t) + ^ fors,t >0 P-a.s., 


(1.4) 


0, m, m- The existence of the deterministic (nonrandom) a.s. limits in ( |1.3[ ) and ( |1.4[ ) follows 
from standard subadditive arguments, which hold in greater generality, for example, when P is an 
arbitrary i.i.d. measure. Define the shape function g by g{s,t) = lim„_>.oo n“^C?([nsJ, [nt J) for 
s,t > 0 P-a.s. The basic properties of g such as concavity and homogeneity are also easily derived 
pni Theorem 2.1]. An interesting matter is the explicit identification of g. Despite much effort, this 
has not been possible except for the exponential and geometric models, which are called exactly 
solvable cases. In this paper, we introduce certain corner growth models with non i.i.d P whose 
shape functions exist and can be determined explicitly. 

We study inhomogeneous generalizations of the exponential and geometric models in which the 
parameters A and q are site-dependent and drawn randomly from an ergodic distribution. More 
specifically, let a = (a„)„gN and b = (bn)neN be stationary random sequences in (0, oo) such that the 
distribution p, of (a, b) is separately ergodic under the map x ti for each fc, Z G N, where tj, is the 
shift (c„)„gN (cn+fc)n€N for k G . In particular, a and b can be independent i.i.d. sequences. 
Suppose that, given (a, b), the weights are conditionally independent and W{i,j) is exponentially 
distributed with rate A = 4- bj for i, j G N. The inhomogeneous geometric model is defined 

similarly except now a and b are sequences in (0,1) and, given (a,b), W{i,j) is geometrically 
distributed with parameter q = aibj for G N. Let us write P for the joint distribution of the 
weights and Pa.b for the joint conditional distribution of the weights given (a, b). To be clear, 
Pa,b is a product measure on and P(il) = JPa,b(S)/i(da,(ib) for any Borel set B c . 
Identifying lF(i, j) with the projection ^ IR+ onto coordinate (i,j). 


Pa.b(fP(b j) ^ x) = e Gi+bj)x Jqj, X ^ 0 


for the exponential model and 


Pa,b(fP(*, j) > fc) = a!lb^ for *, j G N and k e 

for the geometric model. A noteworthy feature that distinguishes these models from the classical 
homogeneous counterparts is the correlations of the weights along the rows and the columns, that 
is, W{i,j) and W{i',j') are not independent under P if j or j = f. 

Our main result is a simple variational description of the shape function. To state it, let a and 
(3 denote the distributions of a = Oi and b = bi, respectively. For the exponential model, 


9{s,t) 


i„t hr 

ze[-g,^] ( Jq 


a{dx) 
X -\- z 



for s, t > 0 , 


(1.5) 


where ij is the left endpoint of the support (the complement of the largest open 77 -null set) of a 
Borel measure 77 on M. (We will also use 77 for the right endpoint of the support). When a and (3 
are uniform measures, we can compute g explicitly. For example, if the supports of a and (3 are 
the interval [1/2, 3/2] then, for s,t > 0, 


, , , t — s ^(t — sV 16st A , / s — t^(s — tV + 16sA _ 

5 (s,t) = slog(^2-f~^- ^ ^ j + f log (^2 + -- )■ (1- 


6 ) 


We obtain similar results for the geometric model in which explicit formulas arise when a and 
P have densities proportional to x 1 —> 1/x. We deduce from (1.51 that g is linear or infinite if 
q = /? = 0, and is differentiable if q -I- /3 > 0. In the latter case, in contrast with (1.3), g may be 
linear close to the axes depending on the behaviors of q and /3 near q and p. More precisely, there 
exist constants 0 ^ ci < C 2 < 00 such that g is strictly concave only inside the cone Ci < s/t < C 2 , 
see Corollary |2.3| This is illustrated through the level set 5 = 1 in Figure |1.1| below. We expect 
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Figure 1.1: The plot oi g = 1 for (1.6) where Ci = 0 and C 2 = QO (left). The plot of g = 1 
and the rays s/t = ci = (—8 + 121og2)/3 « 0.105922 and s/t = C 2 = 4/(9 — 12 log 2) « 
5.863092 for (1.5) when a((ix) = 3a;^dx and/3(dx) = I{is 5 a;s 52 } 4(a;—l)^(ia: (right). 


the finer statistics of G'([ns J, [ntj) as n ^ oo to be qualitatively different in the linear and concave 
sectors. This has been confirmed for large deviation properties in [^. A project currently underway 
is to understand the limit fluctuations. 

A short discussion of the technical aspects of the paper is in order. To calculate the shape 
function, we rely on certain stationary processes with explicit product-form distributions. This 
approach dates back to [27] and is illustrated in |2^ to derive o, which we briefly outline. 
Introduce a parameter z e {q,l/q) and boundary weights {W{k,0),W{0,k) : k e Z+j such that 
: i,j e Z+} are independent, 1F(0, 0) = 0, and W{k,0) and lT(0,fc) are geometrically 
distributed with parameters qjz and gz, respectively. Choose the boundary values in recursion (1.1) 
as G{i, 0) = X;fe=i 0) and G'(0, j) = X;i=i ^(0) k) for z, j e N. Then the resulting model turns 
out to be stationary in the sense that the distributions of the processes {G{i, n) — G{i — l, n) : z e N} 
and {G(n, j) —G(n,/ —1) : j e N} do not depend on n. Consequently, these distributions are product 
measures with geometric marginals. This allows computing the shape function of the stationary 
model and relating it to g via 


5^(1, 1 )= sup max{5^(l-t,0) + g(t, 1),5^(0, l-t)-tg(l,t)} 

t€ [0,1] 


(1.7) 


for z e (g, l/g). Then (1.4) can be extracted from (1.7). The main observation of the present 


work is that, given (a, b), if we choose the parameters of lT(z,0) and VF(0,/) as a^/z and bjZ for 
z G (d, l/j3), we still obtain a stationary model. Then, adapting some arguments from |5S], we also 
arrive at (1.7). To identify g, |29j uses the symmetry g{s,t) = g(t,s), which is not true for the 


inhomogeneous model unless a = (3. For this step, we develop an argument that removes the need 
for symmetry and makes only a few general assumptions on g and gz- 

Literature review. We mention briefly some related results and conjectures beginning with 
the case of i.i.d. P. For the exponential and geometric models. 


g{s,t) = m{s 1) -f 2Vcr^sf, 


( 1 . 8 ) 


wher e m and are the common mean and the variance of the weights. Furthermore, g satisfies 
( jl.Sj ) up to an error of order o{-\/t) as t ], 0 provided that the weights have sufficiently light tail 
|21j . There are also known cases in which g has linear segments. For example, if the weights are 
bounded by 1 and P{W{i,j) = 1) = p > Pc, where Pc is the critical probability for the oriented 
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site percolation, then g{s,t) = s + t in a nontrivial cone [I], [7], [10] ■ Nevertheless, it is expected 
that g is strictly concave and differentiable for a large class of P, for instance, when the weights 
have continuous distributions with enough moments. In our setting, g also enjoys these properties 
under some moment conditions but can be far from (1.8) as (1.6) exemplifies. More recently, 
variational formulas in terms of stationary, integrable cocycles have been developed for g under 
the mild assumption that the weights have finite 2 + e moment for some e > onn. Minimizers of 
these formulas are identified as Busemann functions in m relying on some fixed point results from 
queueing theory m- 

There has also been interest to identify the shape function for non-i.i.d. P. For the exponen¬ 
tial model with columnwise inhomogeneity (i.i.d. a and constant b), |30| obtained a variational 
description of g, which (1.5) includes as a special case. Asymptotics of g near the axes are deter¬ 
mined for more general P in [TSj. Their Theorem 2.4 can be deduced from (1.5). Another direction 
of generalizing the exponential and geometric models is to choose the parameters at site (j,j) as 
A = A(ijn, j jn) and q = Q{ijn, jjn) for some deterministic functions A and Q that encode inhomo¬ 
geneity. Then, under suitable conditions, g can be characterized as the unique monotone viscosity 
solution of a certain Hamilton-Jacobi equation [Ij. 

While we will not take advantage of it in the present paper, we mention that exact solvability 
of the exponential and geometric models goes beyond the explicit limits (1.3) and (1.4). The dis¬ 
tributions of the last-passage times can be expressed as a Fredholm determinant with an explicit 
kernel. Using this, [12] established that G([nsJ, [nt J) has fluctuations of order and converges 
weakly, after suitable rescaling, to the Tracy-Widom GUE distribution. These features are char¬ 
acteristic of the conjectural Kardar-Parisi-Zhang (KPZ) universality class, see survey [5]. More 
generally, Fredholm determinant representations of the distribution of the last-passage times under 
Pa,b have been derived by relating Pa,b to the Schur measures introduced in [23] and, thereby, to 
determinantal point processes with explicit correlation kernels m, [13], m, [IS]. 

Outline. Our results are formally stated in Section [^ We sketch the existence and the basic 
properties of g in Sectio n [3] We discuss stationary versions of the exponential and geometric models 
in Section]^ We prove (1.5) in Section]^ 

Notation and conventions. N = {1,2,3,...} and = {0,1,2,...}. For a; e K, define 
[a:J = max{i e Z : i ^ x} and [x] = min{i e Z : i ^ x}. Also, x v y = inax{x,y} and 
X A y = min{a;,?/} for a:,j/ e K. The imaginary unit is denoted by i. Adjectives increasing and 
decreasing are used in the strict sense. For convenience, we set 1/0 = oo and l/oo = 0. 

Acknowledgement. The author thanks Timo Seppalainen for the problem of computing 
explicit limit shapes and his valuable comments during the preparation of this paper. 


2 Results 

Let E denote the expectation under q. (the distribution of (a, b)). Recall that a = ai and b = bi. 


It is convenient to break (1.5) into the next two theorems. 


Theorem 2.1. Suppose that a + f3 > 0 in the exponential model. Then 


g(s, t) = inf 

ze(-a,§) 


sE 


1 


a + z 


+ tE 


1 


for s,t > 0. 


( 2 . 1 ) 


Hence, g depends on (a, b) only throu gh t he marginal distributions a and f3. Let us write g°‘’^ 
to indicate this. Replacing 2 with —2 in (2.1) reveals that g°‘’^{s,t) = g^’°^{t,s) for s,t > 0, which 


is expected due to the symmetric roles of a and b in the model. In particular, if a and /3 are 
the same then g{s,t) = g(t,s) for s,t > 0. Also, (by dominated convergence) the infimum can be 


taken over [—q,^] in (2.1). When a = § = 0, this interval degenerates to {0} and we expect that 
g{s,t) = sE[l/a] -I- <E[i/ 6] for s,t > 0. Indeed, this is true. 

Theorem 2.2. Suppose that a = (3 = 0 in the exponential model. Then 


gis,t) = sE 


-I- tE 


for s,t > 0. 
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We turn to the concavity and differentiability properties of g. In the case g + /3 > 0, define the 


critical values Ci = 


E[(5+g)-2] 
E[(a-g)-2] 
^-2 


and C 2 = 


E[(a + ^)- 2 ] 


. Note that 0 ^ ci < C 2 < oo. Also, ci = 0 


if and only if E[(o — a) ^] = oo, and C 2 = oo if and only if E[(& — /3) ^] = oo. 


Corollary 2.3. Suppose that a + /3 > 0 in the exponential model. Then 

(a) g{s, t) = s E[(a — g)“^] + t E[(5 + g)~^] for s/t ^ ci. 

(b) g{s,t) = sE[(a + + tF.[{b - for s/t ^ C 2 . 

(c) g{csi + {l-c)s2,cti + {l-c)t2) > cg{si,ti) + { 1 - c)g{s2,t2) force ( 0 , 1 ) and Si,S2,ti,t2 > 0 
such that Cl < si/fi, 52/^2 < C2 and (si, ti) 7^ fc(s2, ^2) for any fc e K. 

(d) g is continuously differentiable. 

By Schwarz inequality, if ci > 0 then E[(a — g)“^] < 00 and if C 2 < 00 then E[(6 — < 00 . 

Hence, g is finite and linear in (s,t) in the regions s/t ^ ci and s/t > C 2 . 


Proof of Corollary \2.3\ Let A{z) = E[(a + z) for 2 ; > —a and B{z) = E[(& — z) for z < p. 
Using dominated convergence, A and B can be differentiated under the expectation. Thus, A'(z) = 
-E[(a + 2)-2], B'{z) =E[(6-z)-2], A"{z) = 2E[(a+z)-3], B"{z) = 2P[{b - z)-% etc. Also, 
define A, B and their derivatives at the endpoints by substituting —a and P for z in the preceding 
formulas. Then, by monotone convergence, the values at the endpoints match the appropriate one¬ 
sided limits, that is, A{—a) = A{z) = E[(a — g)“^], B{p) = lim^|^i3(z) = E[(& — ^)“^], 

and similarly for the derivatives. 

Since A' and B' are increasing and continuous on {—a,p), the derivative z 1 -^ sA'(z) -I- tB'{z) 
is positive if s/t ^ ci, is negative if s/t ^ C 2 and has a unique zero if ci < s/t < C 2 . Hence, (a) 
and (b) follow, and if ci < s/t < C 2 then g{s,t) = sA{z) + tB{z), where z G {—a,p) is the unique 
solution of the equation 


Bfz) _ s 
Afz) ~ t' 


( 2 . 2 ) 


Since —B'/A! is increasing and continuous, it has an increasing inverse (/ defined on (ci,C 2 ). Let 
Si,ti,S 2 ,t 2 be as in (c). Then C(siAi) A C,{s 2 /t 2 ), which implies the strict inequality 


(si + S2)A(z) H- (U + t2)B{z) > SiA(C(si/ti)) -t tiB{(/{si/ti)) 

+ S2A{(/{s2/t2)) A t2B{(/{s2/t2)) (2.3) 

= g{si,tl) + g{s 2 ,t 2 ) 


for any z G {—a,p). Note that ci < (si -I- S2)/(ti + t2) < C2. Setting z = C((si + S2)/(ti + ^2)) 
in ( |2.3| ) yields g{si 4- S 2 U 1 + ^ 2 ) > g{si,ti) + g{s 2 ,t 2 ), and (c) comes from this and homogeneity. 
Since —B'/A! is continuously differentiable with positive derivative (as A",B',B > 0 and A' < 0 


on {—a,P)), by the inverse function theorem, (/ is continuously differentiable as well. Using (2.2), 
we compute the gradient of g for ci < s/t < C 2 as \7g{s,t) = {^{C{s/t)), B{({s/t))), which tends to 
(A(—g), B{—a)) as s/t —> ci and to (A(p),B(p)) as s/t —> C 2 . Hence, (d). □ 


When a and P are uniform distributions, we can compute the infimum in (2.11 explicitly. 


Corollary 2.4. Let A, Z,m > 0. Suppose that a and P are uniform distributions on [A/2, A/2 -I- 1] 
and [A/2, A/2 + m], respectively. Then, for s,t > 0, 


, , s, ( I I It — ms + \/(lt — msY + 4st(A -f Z)(A -I- m) \ 

9(M).pog(^l + - + -- 2.(A,„) - 

t , [ mm ms — It + ^ {It — msY + 4st(A + 1){X + m) 

m ^ 1 A A 2t(A +1) 
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Proof. Since a and /3 are uniform distributions, 


A(z) = E 


1 


a + z 


= y log ( 1 + 


I 


z + Xf2 


B{z) = E 


= — log f 1 
m 


' + A/2 


for z G (—A/2, A/2). We compute the derivatives as 

1 


A'(z) = 


B'{z) = 


(z +A/2)(z +A/2 + Z) (-Z +A/2)(-z +A/2 + m)‘ 

Because = —co and B'{\/2) = co, we have ci = 0 and C 2 = oo. Also, \2.2\ leads to 

(s — t)z^ — (s(A + m) + t(A + l))z + sA(A + 2m)/4 — tA(A + 2Z)/4 = 0. 

It follows from the discriminant formula that the solution in the interval (—A/2, A/2) is 
A s(A + 2?7 i) — t(A + 2/) 

2 s(A + m) + t{X + Z) + -y/(sTO + tiy + 4stA(A + m + 1) 

Inserting this into g{s,t) = sA{z) + tB{z) and some elementary algebra yield the result. 


□ 


The preceding argument can be repeated when Z = 0 or m = 0. In these cases, a and /3 are 
understood as point masses at A/2. For instance, when Z = 0 and m > 0, we obtain 


g{s,t) = 


2sX + ms + {msy + 4stA(A + m) t 


2A(A + TTi) 


+ - log 1 + - + - 

m \ A A 


TO TO ms + ^ (ms)^ + 4stA(A + to) 


2tX 


When Z = 0 and to = 0, we recover (1.3). 

We can also determine g along the diagonal when a and are the same. 

T" 


Corollary 2.5. Suppose that a = (3. Then g{s,s) = 2sE 


for s > 0. 


Proof. We have (a + z) ^ + (a — z) ^ > 2o ^ for |z| ^ a with equality if only if z = 0. Therefore, 


r 1 11 


r 1] 

-^- 

= 2sE 


a + z a — z 


a 


g(s, s) = s inf E 

We only report the analogous results for the geometric model. 
Theorem 2.6. Suppose that aft < 1 in the geometrie model. Then 


□ 


q(s, t) = inf 

z€(a,l/^) 


sE 


af z 

1 — a/z 


+1E 


bz 


1-bz 


for s,t > 0- 


Theorem 2.7. Suppose that a = ft = 1 in the geometric model. Then 


g{s,t) = sE 


1 — a 


+ t E 


1 - b 


for s, t > 0. 


Corollary 2.8. Let q G (0,1) and 0 < Z,to < ^fq. Choose a and ft as the distributions with 
densities proportional to x ^ \/x on the intervals [.y/q — f^/q] and [.y/q — m,.y/q], respectively. 
Then 

s ( IJq I ly-mx + J{ly - mxY + 4a;y(l + to^ - q)(l + l.y/q - q)\ 

g{s, t) = - log 1 + - ---Y'------- 

L \ l — q 1 — ^ 2x{lmy/q — q) I 

t , / mJq m mx — ly^Aj{ly — mxY + Axy{l + rriy/q — q){l + — q)\ 


for s,t > Q, where x = slM, y = tmL, L = log 




Corollary 2.9. Suppose that a = f. Then g{s,s) = 2sE 


and M = log 




■s/g-m 

for s > 0. 
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3 The existence of the shape function 


Lemma 3.1. There exists a deterministic function g : (0,oo)^ —> [0,oo] such that 
liin G([nsJ,[ntJ) ^ for s,t > 0 P-a.s. 

n—^OD n 

Furthermore, g is nondecreasing, homogeneous and concave. 

Here, nondecreasing means that g{s', t') ^ g{s, t) for 0 < s' ^ s and 0 < t' ^ t, and homogeneity 
means that g{cs, ct) = cg{s, t) for s,t,c > 0. In the exponential model, g is finite if a + /3 > 0. This 
is by the standard properties of the stochastic order ED Theorem 1.A3]. Briefly, the i.i.d. measure 
P on under which each W{i,j) is exponentially distributed with rate a + (3 stochastically 
dominates Pa.b and G is a nondecreasing function of the weights. Thus, g{s,t) does not exceed 
the right-hand side of (1.31 with X = a + f3. Similarly, g is finite in the geometric model if aft < 1. 
Extend g to by setting g{0, 0) = 0, g{s, 0) = limt^o g{s, f) and g(0, f) = lims^o g{.s, t) for s, t > 0. 
Lemma |3.1| can be proved using the ergodicity properties of P and superadditivity of the last- 


passage times. As this is quite standard, we will leave out many details. For k,l e Z_|_, let 
be given by 9k,i{uj){i, j) = uj{i + k,j + 1) for i,j e N and w e . Note that P 
is stationary with respect to 6k^i because P(0^J(H)) = EPa.b(^fe = EP,-,^(a),ri(b)(B) = P(B) 

for any Borel set B c . 

Lemma 3.2. P is ergodic with respect to 9k,i for any k,l e N. 

Proof. Suppose 9^\{B) = B for some Borel set B c . For n > 1, let Tn denote the u-algebra 
generated by A„, the collection of W{i,j) with i > k{n — 1) and j > l{n — 1). Then B is in T = 
PlneN Also, T is the tail cr-algebra of the cr-algebras generated by A„ \ A^+i. Because Pa,b is a 
product measure, by Kolmogorov’s 0-1 law, Pa,b(B) e {0,1}. Therefore, P(B) = /r(Pa,b(B) = 1). 
On the other hand, 

(rfe X Tl) {Pa,b(B) = 1} = {Prfca,Tib(B) = 1} = {Pa,h{9kj (B)) = 1} = {Pa,b(B) = 1}. 

Since p is ergodic under Tk x t;, we conclude that P(i?) e {0,1}. □ 


' (iJ)e 7 T 


Proof of Lemma \3.1\ Fix s, t G N and define, for integers 0 ^ m < n 
Z(m, n) = -G((n - m)s, {n - m)t) o 9ms,mt = 


■neHn 


max 

+ l,mt + l,n 


Using the definition and Lemma 3.2 we observe that {Z(m, n) : 0 ^ m < n} is a subadditive process 


that satisfies the hypotheses of Liggett’s subadditive ergodic theorem HZ]. Hence, ^(0,n)/n = 
G{ns, nt)ln converges P-a.s. to a deterministic limit, g{s, t). The existence of the limit for all s, t > 0 
P-a.s. and the claimed properties of g follow as in the case of i.i.d. weights [291 Theorem 2.1]. □ 


4 Stationary distributions of the last-passage increments 

Let us extend the sample space to Now W{i,j) denotes the projection onto coordinate {i,j) 
for i,j G Z_|_. Define the last-passage time G(i,j) through recursion (1.1) but with the boundary 


values G(z,0) = X;fc=i W{k,0) and G(0, j) = X;fe=i W{Q,k) for i,j G N. We then have 
G(m,n) = max > lF(z,B form,nGZ_|_. 

Treno.O.m.n ..“V 

(2j)€7r 


(4.1) 


In the exponential model, for each value of (a, b) such that a„ > a and bn ^ /3 for n G N (which 
holds /r-a.s.) and parameter z G {—a,P), define P^.b the product measure on by 

Plb(bb(i,i) ^x) = exp(-(ai 4- bj)x) Pa b(lT(0,0) = 0) = 1 

Pa.b(l^(b0) >x)= exp(-(a, -t z)x) P^ b(lT(0, j) ^ x) = exp{-{bj - z)x) 
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(4.2) 






for X ^ 0 and i,j e N. When g = /3 = 0, we make definition (4.21 for z = 0. Note that 
the projection of j, onto coordinates is Pa.b- For the geometric model, the construction is 
similar. For 2 e (a, 1//3) and each value of (a, b) such that a„ ^ d and ^ /3 for n G N, the 
measure P^ b is given by 


>k) = a'^b’; 

^ k) = al/z^ 


Plb(^F(0,0) = 0) = l 


k 


Plb(^F(0,j) ^k) = 


(4.3) 


for k G Z+ and i,j G N. When d = /3 = 1, definition (4.3) makes sense for z = 1. 

Introduce the increment variables as I{m,n) = G{m,n) — G{m — l,n) for m > 1 and n ^ 0, 
and J{m,n) = G{m,n) — G{m,n — 1) for m > 0 and n > 1. We capture the stationarity of the 
increments in the following proposition. 

Proposition 4.1. Letk,leZ+. Under 

(a) I{i,l) has the same distribution as W{i,0) forieN. 

(b) J{k,j) has the same distribution as IF(0,j) forj G N. 

(c) The random variables {I{i,l) : i > k} kj {J{k,j) ■ j > 1} are (jointly) independent. 

(1.1) leads to the recursion [521 (2-21)] 

I(rn, n) = I{m, n — 1) — I{m, n—1 )aJ(to — I,n) + 1F (m, n) 

J(rn, n) = J{m — 1, n) — I{n, n — 1) a J{m — 1, n) + lF(m, n) 


(4.4) 


for m,n G N. Proposition 4.1 can be proved via induction using (4.4) and Lemma 4.2 below. We 


will omit the induction argument as it is the same as in 


Lemma 4.2. Let F : 


denote the map {x, y, z) 


H Theorem 2.4]. 

{x — XAy + z,y — xAy + z,XAy). Let 


P be a product measure on K with marginals Pi, P 2 , P 3 . Suppose that one of the following holds. 

(i) Pi, P 2 and P 3 are exponential distributions with rates a, b and a + b, for some a, be (0, 00 ). 

(ii) Pi, P 2 and P 3 are geometric distributions with parameters a, b and ab, for some a,b e (0,1). 
Then P{F~^{B)) = P{B) for any Borel set B czM.^. 


In earlier work |29l Lemma 2.3] and [21 Lemma 4.1], Lemma 4.2 was proved by comparing the 
Laplace transforms of the measures P and P{F~^{-)). We include another proof below. 


Proof of Lemma 14-H We prove (i) only as the proof of (ii) is the discrete version of the same 
argument and is simpler. Observe that F is a bijection on with F~^ = F. It suffices to verify 
the claim for any open set B in By continuity, F~^{B) is also open. Furthermore, F is 
differentiable on the open set {{x,y,z) : x > y or x < y} and its Jacobian equals 1 in absolute 
value. Hence, by the change of variables [221 Theorem 7.26], 


P{F-\B)) = ab{a + b) J 

F-UB) 

= ab{a + b) J 


^—ax—by— 




^—a{x—x Ay+z) — b{y—x Ay+z) — {a+b){x A 


^^dxdydz 




= ab{a + &) J = P{B). □ 


In the exponential and geometric models, respectively, define 


gz{s,t) = sE 
gz{s,t) = sE 


1 


a + z 
a) z 

1 — a/2 


+ tE 


1 


+ tE 


b-z 
bz 


l~bz 


for s, t > 0 and z G [—a, 

for s, t > 0 and z G [d, 1//3]. 





















Lemma 4.3. In the exponential model, let z e (—a,/I) ifa + /3>0, and let z = 0 and assume that 
E[l/a + 1/6] < CO if a = P = 0. In the geometric model, let z e (a, 1//3) if ap < 1, and let z = 1 
and assume that E[a/(1 — a) + 6/(1 — h)\ < co if a = P = 1. Then 

lim ——— = gz(s,t) fors,t^Oin -probability for p-a.e. (a, b). (4-5) 

n—>ao n ' 


In fact, the convergence in (4.5) is Pg^-a-S. for /r-a.e (a, b) provided that a + p > 0 in 
the exponential model and ap < 1 in the geometric model [H Theorem 4.3]. By (1.2), (4.1) 
and nonnegativity of weights, G{m,n) ^ G(rn,n) for m,n e N. Then Lemma 4.3 implies that 
g{s,t) ^ gz{s,t) for any s ,6 0. The main result of this paper is that g{s,t) = iniz gz{s,t). 


Proof of Lemma \4-3\ We will consider the exponential model only, the geometric model is treated 
similarly. Note that G([nsJ, [nt\) = Ji['^s\,j) for s,t ^ 0 and n e N. By 


Proposition 


4.1 


Hence, it suffices to show that 


{J([nsJ,j) : j e N} has the same distribution as {J(0,j) : j e N} under P^.b- 


lim — V /(i, 0) = E 
n—>00 n 


i = l 


1 


a-\- z 


and lim — V J(0,j) = E 

ra —>00 n 


i=i 


in P^ j, for /r-a.s. We will only derive the first limit above, for which we will show that, for z > —a 
and for z = —a when E[(a — a)“^] < oo. 


lim — V /(*, 0) = E 

7,—>nn n 


n^CO Tl 


1 


a + z 


in p-a.s., 


(4.6) 


where is the product measure on the coordinates N x{ 0 } given by Qa(W(j, 0 ) > x) = 

for z e N and x > 0. It suffices to prove the convergence in distribution under /i-a.s. because 

the limit is deterministic. 

The characteristic function of n~^ 2"=! 0) under is given by 


HP 


IX 


n{ai + z) 


-1 


= exp 


£ log (1 


IX 


n{ai + z) 


for X G . 


where the complex logarithm denotes the principal branch. Hence, (4.6) follows if we prove 

1 


lim — > log ( 1 — 

n.—\ 


2=1 


IX 


n{ai + z) 


= ixE 


a + z 


for X G . 


p-a.s., 


Using the bound | log(l + ix)| ^ jxj for x G K and the ergodicity of a, we obtain 


lim sup 

n—»-GO 


2 log ( 1 - 


2=1 


IX 


n{ai + z) 


^ lim — "V - = |x| E 

n—►GO Ti rtz A- 7 


2=1 


Oi + z 


1 


a + z 


for X G . 


p-a.s. 


Therefore, it suffices to prove the following for x G K p-a.s. 

I 


lim — 7 arg ( 1 — 

r?,—►GO \ 


2=1 


IX 


n{ai + z) 


— xE 


a + z 


= lim > arctan 

r?,—j-GO 


2 = 1 


(oi + z) 7 n{ai + z) 


= 0 . 


(4.7) 


Since arctanx = ^ 9(1 + vf) ^du, we can rewrite the second sum above as 
HU" 

TL p 

S I 

2 = 1 


I du X 1 

r du X 1 

r x^x^ dv 

J 1 + n{ai + z) ^ J 

I + n ^ J 

' + x^x^ 
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where we changed the variables via u = vxjn. Pick M > 0. The limsup as n ^ oo of the absolute 
value of the last sum is bounded /r-a.s. by |a:| times 


(ai+z) ^ 


lim 

n—>00 n 


I 


+ x'^v^ 


= E 


(o+z) 


I 


XV 

+ x^v"^ 


where the a.s. convergence is due to the ergodicity of a and the integrability of 


I 


(a+z) t - 2,,2 2 ^ 


+ x'^v'^ a + z 


The last integral is monotone in x^ and vanishes as M ^ co. Hence, (4.71 holds for a; e K /r-a.s. □ 

The next proposition relates to g through a variational formula. 

Proposition 4.4. 


5 z(l,l)= sup max{g^(l - t, 0 ) + g(t, l),g 2 ( 0,1 - t) + g(l,t)} 

te [0,1] 

For z e {—a,P) in the exponential model and for z e (d, l/,0) in the geometric model. 


(4.8) 


Proof. Fix z e {—a,f3) in the exponential model. Since g ^ gz and gz is linear, (4.8) with instead 
of = is immediate. For the opposite inequality, we adapt the argument in [291 Proposition 2.7]. It 
follows from (1.2) and (4.1) that 

G(n, n) = maxmax{G(fc, 0) + G(n — fc + 1, n) o Ok-ip, G(0, k) + G(n, n — k + 1) o 0o,fe-i}- (4.9) 

fce[n] 

Let L e N and consider n > L large enough so that [(i + l)n/L] > [ in/L] for 0 ^ z < L. For any 
k e [n] there exists some 0 ^ i < L such that [ in/L] < k ^ \{i + l)n/L], and the weights are 
nonnegative. Therefore, (|4.9|) implies that 


G(n, n) ^ max max{G([(z + l)n/L], 0) + G([(l — i/L)n\,n) o 9rin/L^ Oj 

0sSi<L ' ’ 

G(0, [(z + l)zz/L]) + G{n, [(1 - i/L)n\) o 0o,[m/Ll}- 
By stationarity of P, we have the following limits in P-probability. 

G([(l-Vi)zzJ,n)o6»fi„/i]^o 

hm -= 5 ( 1 - ^/L, 1 ) 

n—»-00 n 

Ita G(".K1-‘A)"])°»0.I.-/U . ,J_ 1 

n—>00 n 


(4.10) 


(4.11) 


Hence, these limits are P-a.s. and, consequently, Pa,b a.s. g-a.s. if zz —» oo along a suitable sequence 
(nfc)fceN. Also, by Lemma 4.3 there is a subsequence (n'jjkm in N /z-a.s. such that P^b a.s. 


G(f(z + l)n;/Ll,0) 


hm m + 1/2., 0) lim ^ = g^( 0 ,z + l/L) (4.12) 


G(0, [(z + l)</Ll) 


fe —>00 


fc —>00 


rii. 


Because Pa,b is a projection ofP^ j,, we can choose (a,b) such that (4.11) and (4.12) hold P^ ^-a-S. 
Hence, we obtain from (4.10) that 

5 ^( 1 ,1) sS max max{ 52 ((z + 1)/L, 0) + 5(1 - i/L, 1), 5 ^( 0 , (z + 1)/L) + 5 ( 1 ,1 - i/L)} 

0 ^i<L 


^ sup max{g(t, 1 )+ 5 ^(l-t, 0 ), g{l,t) + gz{ 0 ,l - t)} + 


E[(a + z)-i] + E[(&-z)-i] 


Finally, let L ^ 00 . The geometric model is treated similarly. 


□ 
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5 Variational characterization of the shape function 

We now prove Theorems 


2.1 


Theorem 2.2 We begin wit 


and 


2.2 


The assumption a + /3 > 0 is in force until the proof of 


1 computing g on the boundary. Recall that g is extended to the 


boundary of through limits. By homogeneity, it suffices to determine g(l,0) and g(0, 1). 

Lemma 5.1. 


5(1,0) =E 


1 


d + , 


5(0,1) = E 


1 


b + a 


Proof. We have g{l, 0) ^ 52(1,0) = E[(a + 2)“^] for all z e (—a, §). Letting z ] § yields the upper 
bound 5(1,0) ^ E[(a + ^)“^]. Now the lower bound. Let e > 0. By Lemma 3.1 (4.6) and since 
g{bi ^ + e) > 0, there exists (a, b) such that 61 ^ ^ + e and 


G(n, IneJ) , , 

hm -=5(l,e) Pa.b-a.s. 

n—►OO n 


1 


lim — V I(i, 0) = E 
n—>-00 n ^ 


i = l 


1 


(2 + + C 


in -probability. 


(5.1) 

(5.2) 


(Qa is defined immediately after ( |4.6[ )) The distribution of {W{i,\) : 1 ^ i ^ n} under Pa.b 

stochastically dominates the distribution of : 1 ^ i ^ n} under Qa as these distributions 

have product forms and ith marginals are exponentials with rates 0^ + 61 ^ + /3 + e for is \n\. 

Therefore, for cc e K and n > 1/e, 


Pa,b(G(n, [nej) > nx) ^ Pa,b ( lT(i, 1) > na; ] > Qf"’"'' ( ^ /(t,0) > 


nx 


Set X = E[(a + /3 + e) — e and let n —» oo. By (5.1) and (5.2), we obtain 5(1, e) ^ x. Sending 
e ], 0 gives <7(1, 0) 5= E[(a + Computation of 5(0,1) is similar. 


□ 


We now extract g from For this, we will only use the boundary values of g provided in 

Lemma 5.1 and that A(z) = E[(a + z)“^] and B{z) = E[(6—0)“^] are continuous, stricly monotone 
functions on {—a, 13). 

Lemma 5.2. Let r he a positive, continuous function on [0,7r/2]. For z e {—a,/3), 

sup {g{x{ 0 ),y{ 0 )) - g„{x{ 0 ),y{ 0 ))} = 0, 

0^e^7r/2 

where {x{ 0 ),y{ 0 )) = (r(5) cos 0, r(5) sind) for 0^0^ tt/ 2 . 

Proof. We can rewrite ( |4.8| ) as 

A{z) + B{z) = sup {(1 — cot0)A(z) + (7(cot0,1)} V sup {(1 — tan0)B(z) + 5(1,tan0)} 

7r/4sSesS7r/2 0^e^7r/4 


= sup 
71/4^0^71/2 

V sup 
0^0^71/4 


1 ■ 

1 - 


xjd) 

y{0) 

x{ 0 ) 


A{z) + g 


fm 




where we use that x and y are nonzero, respectively, on the intervals [0,7r/4] and [7r/4,7r/2]. Col¬ 
lecting the terms on the right-hand side and using homogeneity, we obtain that 


0 = max 


1 

tv/4:^9^t 7/2 yyft) 


-x{0)A{z) - y{0)B{z) + g{x{0),y{0))}, 


sup —r^ 

0^esS7r/4 ®(u) 


[-x{0)A{z) - y{0)B{z) + g{x{0),y{0))} / 


(5.3) 
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The expressions inside the supremums in (5.3) are continuous functions of 9 over closed intervals. 
Hence, there exists 9^ e [0, 7 r/ 2 ] such that 

0 = -x(9:,)A(z) - y{9;,)B{z) + g{x{9:,),y{9^)) = sup {-x{9)A{z) - y{9)B{z) + g{x{9),y{9))}, 

0sS6)sS7r/2 


where the second equality is due to g ^ g^. 

Corollary 5.3. 


B{z) = snp {-sA{z) + g{s,l)} forze{-a,§). 

0^s<oo 


□ 


(5.4) 


Proof. Let S' > 0. The set {(s, 1) : 0 ^ s ^ S} u {{S,t) : 0 sS t ^ 1 } is the image of a curve 
9 1 -^ {r{9) cos 9, r (9) sin 9) for [ 0 , 7 r/ 2 ] with continuous and positive r. Hence, by Lemma [5. 2 [ 


0 = max{ sup {g(s, 1 ) - ^^(s, 1 )}, sup {g{S,t) - 9z{S,t)}}- 

OsSssSS OsStsSl 


(5.5) 

OsStsSl 

Using homogeneity and Lemma |5.1[ we observe that 

g{S, t) - g^{S, t) = g{S, t) - SA{z) - tB{z) ^ S(g(l, 1/S) - A{z)) -co as S ^ oo. (5.6) 


Hence, the second supremum in (5.5) can be dropped provided that S is sufficiently large, which 
results in 0 = snpQ^g^g{g{s, 1) — (?z(s, 1)}. This equality remains valid if S is replaced with co by 
(5.6) with t = 1. Rearranging terms gives (5.4). □ 

Proof of Theorem \2. 1\ Define the function 7 : M ^ IRu{oo} by 7 ( 5 ) = —g(s, 1) for s > 0 and 
7 ( 5 ) = 00 for s < 0. By Proposition 3.1 7 is nonincreasing, continuous and convex on [0, 00 ) and 
completely determines g. Let 7 * denote the convex conjugate of 7 , that is, 


7 *(x) = supjsa; — 7 ( 5 )} = supjsx — 7 ( 5 )} for x e . 

seE s^O 


(5.7) 


Let / be the function whose graph is the image of the curve z 1 —> {—A{z), B{z)). That is, / is 
defined on the interval (—H(—q), —H(/3)) and is given by the formula f{x) = B o A~"^{—x). By 
Corollary |5.3| 

f{x)= sup {sx — 7 ( 5 )} for X G (—H(—q), —H(^)) (5.8) 


0:$S<OO 


Comparison of (5.7) and (5.8) shows that 7 * coincides with / on (—H(—q), —H(^)). Since 7 is a 
lower semi-continuous, proper convex function on the real line, by the Fenchel-Moreau theorem, 7 
equals the convex conjugate of 7 *, hence. 


7 ( 5 ) = sup{sx — 7 *(x)} for s G . 

ireE 


(5.9) 


To prove the result, we need to show the supremum in (5.9) can be taken over the interval 
(—H(—q), —H(/3)) instead of the real line. It is clear from (5.7) that 7 * is nondecreasing and is 
bounded below by — 7 ( 0 ) = g{0, 1) = B{—a). Since 7 * agrees with / on (—H(—q), —H(/3)), 


H(-q) ^ 7 *(-H(-q)) ^ lim /(x) 

xi—A( — a) 

= lim B o A~^{—x) = lim B{z) = B{—a), 

xi-A{-a) z^-a 


(5.10) 


where we used continuity of A ^ and B. Hence, 7 *(x) = B{—a) for x ^ —A{—a). On the other 
hand, if x > —A{P) = —( 7 ( 1 , 0 ) then 7 *(x) = 00 by (5.7) because 

lim sx — 7(5) = lim s(x -I- g{l, 1/s)) = QO. 

S—>00 S—>00 

Finally, we compute 7 * at —H(/3). Being a convex conjugate, 7 * is lower semi-continuous. Since 
7 * is also nondecreasing, limyi^: 7*(2/) = 'y*{x) for any x G K. Then, proceeding as in (5.10), 

7 *(-H(/ 3 )) = lim /(x) = H(/3). 

xf-A(/3) 
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We conclude that the function x ^ sx — 7 *(x) is increasing for x ^ —A{—a) and is —oo for 
X > —A{j3). Moreover, the left- and right-hand limits agree with the value of the function at 
—A{P) and —A{—a), respectively. Hence, by (5.9|, 


7 (s) = sup {sa; • 

se{-A{-a),-A{§)) 


7 *(a;)} = sup {—sA{z) — B{z)} = — inf {sA{z) + B{z)}, 

ze{-a,§) ^e(-a,l3) 


which implies (2.11. 


□ 


Proof of Theorem \2.^ Introduce d > 0 and let ip : ^ K” denote the map (c„)„gN (cn V(5)„gN. 

Because p commutes with the shift ri, (p{si) and </?(b) are stationary sequences in (0,oo). Moreover, 
for each e N, the distribution ps of ((^(a), (^(b)) is ergodic with respect to x t;. To see 
this, suppose that B = x ti)~^{B) for some /c, Z e N and Borel set B a xK^. Then 
{p X = (v? X X Ti)“i(B)) = (rfc X Ti)-^{{p x Lp)-^{B)). Hence, by the ergodicity 

of p, we get ps{B) = p{{ip x p)-^{B)) e {0,1}. 

Let as and (3s denote the marginal distributions of p{b) and <p(b), respectively. Then as = 
j3_s = S. Applying Theorem 


2.1 


gives 


inf isE 

ze{-S,S) ' 


1 


a V 6 + z 


-f t E 


1 


b V S — z 


(5.11) 


Since Pa,b stochastically dominates P,^(a),¥ 3 (b)) we have for s,t 5= 0. Using 

this and (5.11), we obtain 


5 “’^(s,t)> inf -^sE 
z€(-5.<5) ' 


1 

a V 6' + z 


+ tE 


1 

b V 6' — z 


where we fix d' > S. Because the expression inside the infimum is continuous in z, letting (5 ], 0 
yields > sE[(a v + iE[(& v foi' s,t > 0. Then, by monotone convergence, 

letting i5' —> 0 results in 


g^’P{sA)^sP 


-t t E 


for s, t > 0. 


The opposite inequality is noted after Lemma |4.3| 


□ 
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